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ABSTRACT 

We present observations of interstellar rubidium toward o Per, £ Per, AE 
Aur, HD 147889, x Oph, ( Oph, and 20 Aql. Theory suggests that stable 85 Rb 
and long-lived 87 Rb are produced predominantly by high-mass stars, through a 
combination of the weak s- and r-processes. The 85 Rb/ 87 Rb ratio was determined 
from measurements of the Rb I line at 7800 A and was compared to the solar 
system meteoritic ratio of 2.59. Within 1-cr uncertainties all directions except 
HD 147889 have Rb isotope ratios consistent with the solar system value. The 
ratio toward HD 147889 is much lower than the meteoritic value and similar to 



that toward p Oph A (jFederman et al.l 12004 ); both lines of sight probe the Rho 



Ophiuchus Molecular Cloud. The earlier result was attributed to a deficit of r- 
processed 85 Rb. Our larger sample suggests instead that 87 Rb is enhanced in these 
two lines of sight. When the total elemental abundance of Rb is compared to the 
K elemental abundance, the interstellar Rb/K ratio is significantly lower than the 
meteoritic ratio for all the sight lines in this study. Available interstellar samples 
for other s- and r- process elements are used to help interpret these results. 

Subject headings: ISM: abundances — ISM: atoms — stars: individual (o Persei, 
C Persei, AE Aurigae, HD 147889, x Ophiuchi, ( Ophiuchi, 20 Aquilae) 



1. Introduction 

Rubidium isotopes 85 Rb and 87 Rb are products of the neutron capture slow (s-) and rapid 
[r-) processes. ( 87 Rb is unstable with the long halflife of 4.75 x 10 10 years.) The s-process 
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contributions come from two types of stars: the so-called main s-process operates in low 
mass asymptotic giant branch (AGB) stars with their products ejected into the interstellar 
medium (ISM) by the stellar wind; and the weak s-process occurs in the He- and C-burning 
shells of massive stars and distributed into the ISM primarily, it is thought, by a star's 
terminal supernova explosion. The r-process most likely occurs at the time of the supernova 
explosion. Of primary relevance to studies of the Rb elemental and isotopic abundances 
beyond the solar system is the division of responsbility for Rb synthesis between low-mass 
long-lived and high-mass short-lived stars or, equivalently, the relative contributions from 
the main s-process and the combination of the weak s-process and the r-process. 

Dissection of the solar system abundances for Rb and adjacent elements enables frac- 
tional contributions of the weak and main s-processes to be estimated. Then, the r-process 
contribution is obtained as the differences between the solar system abundance and the sum 
of the weak and main s-pr o cess contributions. The most recent exercise of this kind is 
that reported by iHeil et al.l (120081 ) who estimated the fractional contributions of the main 
s-process to the solar system abundances of 85 Rb and 87 Rb to be 0.17 and 0.24, respectively. 
In other words, the relative contributions of low mass and high mass stars is 0.17/0.83 = 0.20 
for 85 Rb and 0.24/0.76 = 0.32 for 87 Rb. [Heil et all (120081 ) give the weak s-process fractional 
contributions as 0.24 for 85 Rb and 0.46 for 87 Rb and, then by subtraction of the two s-process 
contributions, the r-process contributions are 0.59 for 85 Rb and 0.30 for 87 Rb.] Earlier dis- 
sections of the solar system abundances have g iven somewhat d i fferen t results for the relative 
contributions of low and high mass stars, e.g.. lArlandini et al.l (Il999l ) gave 0.16/0.84 = 0.19 
for 85 Rb and 0.35/0.65 = 0.54 for 87 Rb indicating a more pronou nced role of the massive 
stars in the synthesis of 87 Rb (relative to 85 Rb) than suggested by IHeil et al.l (120081 ) . Most 
probably, higher weight should be given to Heil et al.'s result because it was based in part 
on new and accurate measurements of neutron capture cross sections, especially for the two 
Rb isotopes. Our study of interstellar Rb was undertaken to set observational constraints 
on the proposed mechanisms of Rb nucleosynthesis. 



Rubidium has been analyzed in both une volved and evolved stars. lLambert &: Luck 



( 119761 ) sought the Rb isotope ratio in Arcturus. iGratton fc Snedenl (119941) studied the abun- 
dances of neutron-rich elements in metal-poor stars, while iTomkin fc Lambertl (Il999l ) con- 
sidered rubidium in metal-deficient disk and halo stars. Others have examined rubidium 



abundances in AGB stars (lLambert et al.lll995t lAbia fc Wallersteinlll998l ; lAbia et al.ll200ll ). 
However, rubidium is difficult to measure in stellar atmospheres. The Rb I lines at 7800 
and 7947 A are blended with nearby lines and stellar broadening prevents the hyperfine 
and isotopic components from being resolved. Intrinsically narrow interstellar lines are not 
hindered by these complications. 
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The detection of rubidium in the ISM was first reported by lJura fc Smith! (Il98ll ) from 
absorption toward the bright star ( Oph. They used the Rb I 5s 2 Si/ 2 -5p 2 P%/2 line at 
7800.29 A. This is the strongest resonance line of neutral rubidium and it is used for our 
study as well, even though most rubidium in diffuse interstellar clouds is ionized because it 
has a low ionization potential of 4.177 eV. Rubidium's overall low abundance requires the 
use of its strongest line. 



The early detection by lJura fc Smith! (jl98ll ) could not be confir med by IFederman et al. 

( 1985 ). who also could not detect Rb I toward o Per nor £ Per. Federman et al. ( 19851 ) 
calculated a typical 3-a upper limit in column density of ~ 5 x 10 9 cm -2 for these lines of 



sight. Interst e 



lar a bsorption toward these three stars is revisited in this study. More recently, 



Gredel et al.l (120011 ) provided the first firm detections of interstellar rubidium toward Cyg 
OB2 No. 12 and No. 5. 

The solar system isotope ratio is found f rom carbonaceou s chon drite meteorites, with 
a ratio of 85 Rb/ 87 Rb = 2.59 (ILoddersI 120031 ) . Federman et al.l (120041 ) found an interstellar 
isotope ratio toward p Oph A of 1.21, which differs significantly from the meteoritic value. 
This was the first relatively precise determination of the rubidium isotope ratio for extrasolar 
gas. In a recent paper, Kawanomoto et al. (2009) studied the line of sight toward HD 169454 
and derived an isotope ratio that is consistent with t he solar system value. The present study 
expands upon the results of IFederman et al.l (120041 ) and Kawanomoto et al. (2009) on the 
stable 85 Rb and long-lived 87 Rb isotopes. 

We sought Rb absorption along seven lines of sight toward o Per, £ Per, AE Aur, 
HD 147889, x Oph, ( Oph, and 20 Aql. These lines of sight were selected because they 
sample different directions in the solar neighborhood, and because other relevant data noted 
below are available for them. Furthermore, the cloud structure is simple with only one or 
two apparent components per line of sight. We explore 85 Rb/ 87 Rb ratios and compare the 
results to those for p Oph A. We also compare the elemental abundance ratio of Rb/K, as 
well as the 85 Rb/K and 87 Rb/K ratios, to examine whether the 85 Rb or 87 Rb abundances 
are lower or higher relative to the abundances seen in meteorites. Section 2 summarizes the 
high-resolution observations and data reduction techniques on the Rb I A7800 line. Section 
3 presents the analysis of the data and Section 4 the results, including confidence limits on 
the column densities. Section 5 uses the results to discuss the implications for interstellar 
rubidium and its nucleosynthetic production site(s). Finally, Section 6 gives a summary and 
conclusion. 
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2. Observations and Data Reduction 

We observed the star s with the Robert G. Tull Spectrograph - formerly 2dcoude spectro- 



graph - (iTull et al.lll995l ) on the Harlan J. Smith 2.7 m telescope at McDonald Observatory 
from 2005 to 2007. We used the high-resolution mode with echelle grating El to acquire 
spectra on a 2048 x 2048 Tektronix CCD. Depending on the central wavelength, there were 
either 12 or 13 orders imaged onto the CCD and each order was around 30 A wide. The grat- 
ing was centered near 7298 A to capture the Rb I A7800 feature in either the eighth or ninth 
order. A 145 /um slit gave a resolving power of 175,000 (2.9 pixels per resolution element), as 
determined from the width of lines in the Th-Ar comparison spectra. The resolving power 
varied from 172,000 to 184,000 from night to night and from epoch to epoch. For calibration, 
we obtained bias and flat-field exposures at the beginning of each night. Furthermore, we 
acquired a set of 30-minute long dark current exposures on the first night of the observing 
runs. The dark images were used to determine the rate that thermal electrons accumulated 
on the CCD. Every 2-3 hours throughout the night we obtained a Th-Ar comparison spec- 
trum for wavelength calibration. Also each night we observed an unreddened star such as 
a Leo or 7 Cas centered on the slit to see if there were any CCD blemishes not removed 
by the flat-fielding process. There was a seven-pixel wide CCD detector glitch in the 2005 
December run that affected the spectra of £ Per. We were unable to remove the glitch and 
subsequently those spectra were not used. In all other runs the glitch did not fall on the Rb I 
feature. This was due to the fact that we shifted the center position of the detector slightly 
each night to minimize the chance of the feature falling on corrupted pixels. Individual ex- 
posures were limited to 30 m to minimize the presence of cosmic rays. Table 1 lists the stellar 
and observational parameters: HD number, name, spectral type, V magnitude, and distance 
in pc (obtained from the SIMBAD database, operated at CDS, Strasbourg, France), dates 
observed, exposure time, and unreddened stars observed. 

The data were reduced using standard routines within the IRAF environment. Master 
bias and dark images were subtracted from all the raw object, comparison, and flat images. 
These images were processed further to apply the overscan strip correction. Cosmic rays were 
eliminated in the object and comparison images and these images, along with the master 
flat, were fit with a low-order polynomial that subtracted the scattered light from each order. 
The master flat was normalized to unity and divided into each object and comparison image 
and the light from each aperture was summed to form one-dimensional spectra. At this 
point, the spectra of the program and unreddened stars were compared to check for CCD 
glitches or artifacts. Typically three or more Th-Ar lines in the comparison spectra were 
identified per order for the wavelength solution. The rms deviation averaged 10 -4 A by 
fitting the x and y directions to second order polynomials. Over the entire spectrum, there 
were few moderately strong lines. Two comparison spectra were assigned to each spectrum 
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and it was corrected for dispersion. Each spectrum was then Doppler-corrected and all were 
summed to obtain a spectrum for each night. Each night's summed spectra were compared 
to make sure there were no discrepancies in the width or depth of the Rb I feature. None 
were present. Finally, the individual exposures for each line of sight were summed to produce 
seven spectra. 

Table 2 gives the signal-to-noise ratio per pixel and the equivalent width ( W\) of the 
line divided by its uncertainty (a). The equivalent width was measured using the deblending 
function in SPLOT of IRAF, where Voigt profiles were fit to the features. The uncertainty 
was calculated by dividing the full width at half maximum (FWHM) of the Rb I feature by 
the average signal-to- noise (S/N) in the stellar continuum. The goal was to achieve values 
for W\/a of about 10 so that the Rb I profiles could be fit with confidence. 



3. Analysis 



The Rb I spectra were synthesized using a C++ program (J. Zsargo, modified by 
Knauth, hereafter called RbFits) to fit the profile, taking the instrumental width into ac- 
count. RbFits took as input Vlsr, the 6-value, and the column density, N. In order to 
determine the starting input velocities and 6-values for the RbFits program, K I parameters 
were used, as K I is expected to trace Rb I in diffuse interstellar gas. Whenever possible, 
results from the weak K I lines at 4044 and 4047 A, whose values of W\ are similar to 
those of Rb I, were adopted. Occasionally the strong line at A7699 was used to find the 
component structure, but in all cases the weak lines provided K I column densities. It is also 
worth noting that the K I lines have unresolved hyperfine components and are not affected 
by the presence of isotopes; thus their profiles ar e relatively simple. In fact, all the alkalis 



are likely to coexist in diffuse interstellar clouds (I Welty & Hobbsl l200ll ; iKnauth et al.l 12003 



Federman et al.l 120041 ); if the K I did not trace Rb I particularly well, 7 Li I and 6 Li I were 
used as a guide. The rubidium absorption toward \ Oph and 20 Aql incorporated starting 
values from the Li I data; all other lines of sight began their fits with K I velocities and 
6-values. The velocities reported below are based on fits to the Rb I profile. These fina l 
Rb I velocities agree well with those determined for Li I and K I (e.g., IKnauth et al.ll2003l ). 
The RbFits program fit Voigt profiles to the features taking the hy perfine structu re of both 
85 Rb I and 87 Rb I into account. The wavelengths and /-values from iMortonl (120001 ) are given 
in Table 3. 

When fitting the profile, RbFits calculated the lowest x 2 by modifying the input values. 
Also, each input parameter of each cloud was either allowed to be held fixed or to vary. The 
goal was to obtain a synthesized profile where the velocity and 6-values for both isotopes 
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were equal, and the column densities of each cloud formed an isotope ratio at least similar to 
the solar system value, 2.59. The inferred isotope ratio did not depend on the initial value 
of 2.0 or 2.5. At the same time, a check was performed on the reduced x 2 1 which ideally 
should be ~1.0. For lines of sight with one cloud, there were 43 degrees of freedom; there 
were 51 degrees of freedom for lines of sight with two clouds (one degree per data point of 
the feature minus the number of parameters being fit). 

The simplest lines of sight were those that had only one clou d component, ( Per , 
HD147889, and x Oph. The K I template for C Pe r was taken from iKnauth et al.l (120001 ). 
and that for x Oph was from IKnauth et al.l (120031 ). An HD 147889 K I spectrum for the 
line at 7699 A was obtained from our unpublished results and was fitted using a FORTRAN 
InterStellar MODelling program, ISMOD (Sheffer, unpublished), because there was no cloud 
structure available in the literature. Using velocities, 6-values, and equivalent widths esti- 
mated from IRAF routines, ISMOD fit a profile to the HD 147889 data that included unre- 
solved hyperfine structure. The results from the main component were used as the starting 
point for RbFits. 



Federman et al.l (120041 ) performed a synthesis with RbFits on p Oph A. The first step 



in using the RbFits program was to reproduce the results from that study. First, all six 
parameters were allowed to vary. From this output, both velocities were set to 1.96 km 
s _1 and both 6-values were set to 1.43 km s _1 from the solution of 87 Rb. RbFits was run 
again with these four values held fixed, only allowing the colum n density to var y. The 



resulting solution had a reduced x 2 °f 1-57, similar to the 1.43 of iFederman et al.l (120041 ) . 
The 87 Rb/ 85 Rb isotope rati o for p Oph A was com puted to be 1.22 ± 0.30, consistent with 
the 1.21 ± 0.30 obtained by lFederman et al.l (120041 ) . Furthermore, a plot of the spectrum of 
p Oph A was created with the data and fit. A visual inspection was performed to make sure 
that the residuals (data minus fit) were the same inside and outside the profile. 

The above technique was used to obtain syntheses for the seven new lines of sight. The 
basic routine was to start by varying everything and then fixing more and more parameters 
until the solution had a reduced x 2 around 1.0, the residuals were consistent along the profile, 
and the 87 Rb/ 85 Rb isotope ratio was reasonable. 



3.1. One Component Fits 



The start i ng in put values for x Oph were taken from the K I and 7 Li I values of 



Knauth et al.l (120031 ). The input velocity and b- value were 0.20 and 1.00 km s 1 , respec- 



tively. Because the lines are so weak, the input column density was derived from the ratio 
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of equivalent width to column density for p Oph A using the simple equation 



W x (p Oph A) _ W x {x Oph) 
N(p Oph A) N( X Oph) ' 



where W\(X) is the equivalent width of the Rb I line and N is the total column density of 
87 Rb and 85 Rb combined. Once the total column density was determined, it was divided so 
that the isotope ratio would be similar to that of the solar system. Thus, 67% of N(x Oph) 
was 85 Rb and 33% of it was 87 Rb, using a ratio of 2:1 that is midway be tween the results 



for p Oph A (IFederman et al.l 12004 ) and the solar system (lLoddersll2003l ). These were the 
column densities used for the starting point in RbFits. x Oph has the highest reduced x 2 i 
but the synthesized profile fits the data well. This "ratio" technique was also used for HD 
147889. 

The profile for ( Per was synt hesized in thr e e step s, following the above procedure. 
The starting values were taken from lKnauth et al.l (120001 ). During the second step, however, 
instead of using the velocit y output from 87 Rb I, the velocity of 85 Rb I was more similar to 
that of K I and 7 Li I from iKnauth et al.l fl2000[ ). After using the velocity and 6-value from 
the solution of 85 Rb I, the rest of the synthesis was performed as described above. Figure 1 
displays the fitted spectra toward ( Per, x Oph, and HD 147889. 



3.2. Two Component Fits 



The profiles of the four lines of sight with two components were more difficult to syn- 
thesize. The starting values did not affect the outcome of the one-component fits, whereas 
the simple "ratio" technique of dividing equivalent widths by column densities was no longer 
valid because the profiles were more complex, each having eight components (two hyperfine 
components for each isotope for each cloud). The program needed better starting values 
from other data (velocity separations, etc.), especially for the column densities. There are 
two ways to solve this problem. 

For the first way, a list of abundances of Rb I and K I was created from lines of sight 
with one componen t: x Oph, C P er > HD 147889, p Oph A, and the two lines of sight from 
Gredel et al.l (12001 ) Cvg OB No. 5 and No. 12 . The K I colu mn densities were obtained 
fj200oh (C Perl. lKnauth et all (12003^ (x Oph). IFederman et al.l koo4 ) (p 



from Knauth et al 



Oph A, Cyg OB No. 5 and No. 12), and derived from unpublished observations (HD 147889). 
An average of the iV(Rb T)/N(K I) ratios for each line of sight yielded 9.5 x 10~ 4 . At this 
preliminary stage of the analysis for AE Aur, the component structure and iV(K I) come from 
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the synthesis of A7699 from unpublished data. For all other directions, N(K I) is obtained 
from the weak line at A4044. Using the above average and knowing K I for the four lines of 
sight in the current study, one could work backwards to find iV(Rb I). Then usi ng ,/V(Rb I) as 



the to tal column density for a given line of sight and the cloud fractions from iKnauth et al. 



(120031 ) for o Per, £ Oph, and 20 Aql as well as our results for AE Aur, the column densities 
for each cloud were obtained. 



to 



For example, o Per h as a total K I column density of 1.04 x 10 cm . Then, according 



Federman et al.l (120041 ) and assuming equal depleti on of Rb and K although the elements 



have somewhat different condensation temperatures (ILoddersI 120031 ). 



A g (Rb) 

L A,(K) J 



iV(Rb I) 



_ N(K I) J |_ G(K I) 



G(Rb I) 



q(K I) 
a(Rb I)_ 



(2) 



where A g (IVo)/ A g (K) is the elemental abundance ratio, G(X) is the photoionization rate 
corrected for grain attenuation, and a(X) is the rate coefficien t for radiative recombi na- 
tion. We adopted the same sources for the atomic data used by iFederman et al.l (120041 ) to 



infer photoionization rates for Rb I of 3 .42 x 10 12 s 1 and for K I of 8.67 x 10 12 s 1 . 



For Rb I, the theoretical calculations of Weis 



reit 



( 1972) giv e cross sections for 1150-1250 



A, which were scaled to the measurements of 
The cross sections for K I are taken from the measurements of 



Marr &: Creek! (1 19681) at longer wave l engths 



Hudson fc Carter! (Il96fih 



Hudson fc Carter! (119671 ) , iMarr fc Creek! (119681 ) , and ISandner et al.l (119811 ) . For wavelengths 
below 1150 A, there is no differential attenuation because the ionization potentials are sim- 
ilar (4.18 eV for Rb I and 4.34 eV for K I). Based on these cross-sections, a 13% correction 
was applied to G(Rb I) for A < 1150 A. The ra te coefficients for radiative electron recom 



bination are similar, w ithin 10% of one another (!Wandll985l ; iPequignot fc Aldrovandilll986 



Wane fc Aymarl 119871 ) , and so the ratio of the rate coefficients is 
photoionization rates is 0.394. 



1. The ratio of the 



For o Per we can solve for iV(Rb I), which is 2.51 x 10 9 cm" 2 . Then, we mul tiplied 
this column density by the percentage in each cloud [taken from lKnauth et al.l (120031 )] . 14% 
and 86% for the clouds at 3.96 and 6.98 km s -1 , respectively. Assuming the ratio in each 
component would be similar to the solar system value, the relative amounts of the isotopes 
were separated. The result was that four starting column densities were created. 

The second way to obtain starting column densities was to perform a one component fit 
on the data, as in the case of ( Oph. This one component fit provided a coarse estimate of a 
weighted velocity, b- value, and column density for each isotope. The isotopes were then split 
into their resultant components according to IKnauth et al.l ( 20031 ). The resulting spectra 
with two components are shown in Figure 2. 
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4. Results 

4.1. The Rubidium Isotope Ratio 

The Rb isotope ratio is expected to be independent of the two principal effects that 
hamper the extraction of the Rb/H ratio from the Rb I equivalent widths, namely, the large 
degree of ionization of Rb to Rb + , and the loss of Rb onto grains. Table 4 shows the final 
solutions for the lines of sight in our survey. It lists the resulting velocities, column densities, 
b- values, 85 Rb/ 87 Rb ratios, and reduced % 2 s, calculated from the RbFits program. The b- 
values range from 0.45 to 1.46 km s -1 ; while the larger values may indicate the presence of 
unresolved structure, similar values have been reported previously (e.g., Knauth et al. 2003). 
The reduced x 2 s are & U around 1.0. The uncertainties for each cloud were calculated using 
the FWHM, derived from the b- value, and S/N of each spectrum. For such weak lines, the 
same fractional uncertainty applies to the column densities also. These were propagated in 
the usual way to place uncertainties on the inferred isotope ratios. When 2-a was greater 
than the equivalent widths for the 87 Rb line, only lower limits to the ratio are possible. All 
lines of sight are consistent with the solar system value, 2.59, as determined from meteorites, 
except for HD 147889 and the previously measured p Oph A (Federman et al. 2004). 

We investigated the sensitivity of our results for HD 147889 by attempting to fit the 
Rb I profile with the solar system ratio. Various combinations of keeping Vlsr and 6-value 
fixed were attempted. All gave similar results, with reduced x 2 s from 1.62 to 1.81 compared 
to our preferred synthesis with 1.44. We performed analyses of the Rb isotope ratio toward 
HD 147889 of a single component fit in two additional ways. In the first analysis all profile 
synthesis parameters were allowed to vary. For the second analysis only the 85 Rb and 87 Rb 
column densities were allowed to vary. For the second synthesis, the velocity and the Doppler 
broadening parameters were fixed to those determined from the profile fit to the measured 
interstellar K I line. This led to an improved reduced x 2 over that of the freely varying 
profile synthesis. Application of an F-test (Lupton 1993) shows the improvement in x 2 is 
justified to the 68% confidence level. In order to test whether our single component fit of the 
Rb isotope ratio is potentially two components, a subsequent analysis was performed. This 
analysis led to an overall improvement in the reduced y 2 . An F-test shows that a second 
component is justified at the 80% confidence level. The synthesis based on changing Vlsr 
appears in the lower right panel of Figure 1. The fit on the blue side is poorer and below 
the data points, which is evidenced by the residuals within the left side of the profile not 
appearing random. 
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4.2. Comparison with Elemental Potassium Abundance 



Neutral Rb is expected to follow the spatial trend of neutral K in interstellar clouds be- 
cause they have similar ionization potentials and chemical properties, and so a determination 
of the interstellar elemental Rb/K ratio from Eqn. 2 compared to the meteoritic value may 
provide new insights into the production of n eutron capture elem ents. High-resolution sur - 
veys hav e been performed for interstellar Li I (IKnauth et al.ll2003l) . Na I (jWelty et al.lll994j ). 
and K I ( 1 Welty fc Hobbsll200ll ) absorption. IKnauth et al.l (120031 ) conc lude that the elementa l 
Li/K abundance ratio is consistent with the solar system value, and lWelty fc Hobbd (120011 ) 
also confirm the essentially linear relationship between the total column densities of Li I, 
Na I, and K I. It seems that Li, Na, and K are likely depleted by compar able amounts in 
the r elatively cool, dense regions where the neutral species are concentrated (jWelty fc Hobbs 
20011 ). Thus we apply Eqn. 2 with confidence. 



Using this equation Rb/K, as well as 85 Rb/K and 87 R b/K, were determined and the 
results are presente d in Table 5. The results for p Oph A (IFederman et al.l 120041 ) and the 
solar system values (lLoddersll2003l ) are shown for comparison. The K I column densities were 
obtained from the weak line at 4044 A, except for AE Aur, where we adopted the unpublishe d 
A4047 results of P. Boisse. The K I value for HD 147889 comes from lWeltv fc Hobbsl (l200lh . 



The seven new lines of sight given in Table 5 all show a deficit in the ratios relative to 
potassium, with the possible exception of 87 Rb/K toward p Oph A. Relative to the meteoritic 
ratios, we find Rb/K=34%, 85 Rb/K=32%, and 87 Rb/K=36%. These ratios are inferred from 
Eqn. (2), which does not allow for different depletions onto grains. As mentioned below, the 
condensation temperature for Rb is lower than that for K; if only depletion mattered, one 
would expect a higher Rb/K ratio, contrary to the observational results. Note that 85 Rb/K 
toward o Per seems to be larger than the other ratios, and that the anomalously large ratios 
of 87 Rb/K (HD 147889 and p Oph A) are not included in the above averages. 



Discussion 



5.1. Comparison with Previous Results 



The first limits on interstellar iV(Rb I) come from IFederman et al.l (119851 ). who searched 
for absorption toward the stars o Per, ( Per, and C Oph, but did not detect any to a 3-cr 
limit of < 1.5 mA. The resulting l imits on 7V(Rb I) and Rb/K are found to be typically < 



5 x 10 9 cm 2 and < 1.5 x 10 3 (IFederman et al.l 120041 ) . In all cases, our detections yield 
values 30 to 50% of the upper limits. 
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Gredel et al.l (120011 ) found column densities toward Cyg OB2 No. 5 and No. 12 of 



N(Rb I) = (7 ± 2) x 10 9 cm" 2 and N(Rb I) = (13 ± 2) x 10 9 cm~ 2 , respectively . Com - 



bining these results w ith the high- resolution K I A7699 spectra of iMcCall et al.l (120021 ). 
Federman et al.l (120041 ) derived values for N(K I) of the three main molecular components to 
be 9.4(7.7) x 10 11 , 8.1(12) x 10 11 , and 6.9(16) x 10 11 cm" 2 for the gas toward Cyg OB2 No. 
5(12). Thus, Rb/K is 14 x 10~ 4 for No. 5 and 12 x 10~ 4 for No. 12. Rb/K for p Oph A 
is (13 ± 3) x 10 _4 (IFederman et al.ll2004l ). These three lines of sight appear to have Rb/K 
values about a factor of two higher than the average of our eleven clouds, (6.1 ±1.2) x 
10 -4 . However, significant uncertainty remains in extracting elemental abundances toward 
the stars in Cyg OB2 because a blended, optically th i ck lin e was used to obtain the K I 
column density. Our analyses, and those of ISnow et al.l (120081 ). for HD 147889 indicate that 
the weaker line at 4044 A yields a higher K I column density, leading to a smaller Rb/K 
ratio. Furthermore, the difference for the gas toward p Oph A is within ~2-a of the results 
for our other lines of sight. Thus, all interstellar results appear smaller than the solar system 
value of (18 ± 4) x 10" 4 . 



From the Rb I detection toward p Oph A, IFederman et al.l (120041 ) suggested that there 
was a deficit of r-processed material because 87 Rb/K seemed to have the meteoritic value, 
while 85 Rb/K was considerably less than the solar system value. This is also true toward 
HD 147889. These two lines of sight are near one another; the anomalously low 85 Rb/ 87 Rb 
ratios seem to arise from sampling the Rho Ophiuchus Molecular Cloud. For the other lines 
of sight, our sample indicates that both 85 Rb and 87 Rb have a deficit relative to K in the ISM 
of ~66%. However, the gas in the Rho Ophiuchus Molecular Cloud appears to be enriched 
in 87 Rb. 

Kawanomoto et al. (2009) recently published results for the direction toward HD 169454. 
Their isotope ratio ( 85 Rb/ 87 Rb > 2.4) is consistent with the meteoritic value, which again 
highlights the unusual values that we found for the Rho Ophiuchus Molecular Cloud. Their 
Rb/K ratio also reveals a deficit compared to the meteoritic ratio. It is worth noting that 
Kawanomoto et al. make comparisons with a solar system ratio at its formation by account- 
ing for the decay of 87 Rb over 4.55 Gyr. Their estimate of the original ratio is 2.43. 



5.2. Interstellar Results on Neutron-Capture Elements 

The elemental abundances of Ge, Kr, Cd, and Sn seen in our sample of diffuse clouds, 
for which sufficient data exist, may provide an answer, because they are also synthesized 
through neutron capture. For each of these elements, the dominant ion is observed. In the 
solar system about 10% of Ge, 20% of Kr, 50% of Cd, and 70% of Sn are produced via the 
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main s-process (IBeer et al.lll992l ; iRaiteri et al.lll993l ; lArlandini et al.lll999l ; iHeil et al.ll2008l ); 

massive stars through a combination of weak s-process and r-pro cess are respons i ble fo r 
the remaining fractions. When co rrected for the /-value adopted by ICartledge et al.l (120061 ) , 
Cardelli. Savage, fc Ebbetl Jl99lh obtain N(Ge)/N(R) of (4. 9 ± 0.5) x 1Q- 10 for C Oph, while 



Loddersl ( 2003 ) gives (4.2 ± 0.5) x 10 9 for the solar system. Cartledge. Meyer, fc Lauroesch 



(120031 ) find JV(Kr)/iV(H) to be (9.33 ± 1.14) x 10~ 10 for C P er and (7 . 94 ±0 .97) x 10~ 10 for C 
Oph, whereas the theoretical s olar system abund ance from Lodders ( 20031 ) is (19.05 ± 3.81) 
x 10- 10 . Sn was measured bv ISofia et all Jl999h . who found iV(Sn)/7V(H) of (7.41 ± 3.18) 
x 10~ n for C Per, (8.13 ± 0.98) x 10" 11 for x Oph, (9.12 ± 2.65) x 10~ n for C Oph, and 
(5.75 ± 95) x 10~ n for p Oph A. The solar system abundance is (12.88 ± 1.24) x 10" 11 
( jLoddersll2003l ). Finally, there is one line of sight amo ng our sample tha t has measured Cd, 
C Oph, where iV(Cd)/7V(H) is (5.50 ± 1.60) x l O" 11 dSofia et ali ll999h. The solar system 
value for this element is (5.50 ± 0.38) x 10~ n JLoddersI bool T Taking into account the 
uncertainties, all these values are likely to be lower than the solar system values. 

For these relatively dense lines of sight, depletion onto grains, not nucleosynthetic pro- 
cesses, could be the cause for the deficits. It is better to consider the gas phase abundances 
for low density, warm gas. The condensatio n temperature s for Ge, Kr, Cd, Sn, and Rb are 



734, 5 2, 652, 704, and 800 K, respectively (ILoddersI 120031 ). The survey bv ICartledge et al. 



( 120061 ) indicates an interstellar Ge abundance that is 25% the meteoritic value. With such a 
low condensation temperature, Kr shoul d not be depleted toward any line of sig ht. However, 
the average abundance of Kr found by ICartledge. Mever. fe Lauroeschl (120031 ) is ~50% of 
the solar system value. The results on Cd and Sn (ISofia et al.lll999l ) reveal a different trend. 
For low density lines of sight, neither element shows any depletion onto grains. There is 
even a hint that the Sn abundance for these directions is greater than the solar system value. 
The deficit of Rb relative to K cannot be explained by depletion onto grains resulting from 
a higher condensation temper ature, be c ause Li, Na, and K have even higher values (1135, 
953, and 1001 K) according to Loddersl (2003), and the Li/K and Na/K elemental ratios are 
similar to meteoritic values ( Welty fc Hobbsl 2001 ; Knauth et al. 2003 ). 



The depletion patterns of Ge, Kr, and Rb on the one hand, and Cd and Sn on th e 
other, cannot be attributed to imprecise oscillator strengths. According to IMortonl (120031 ). 
the /-values used to study interstellar Kr I, Rb I, and Cd II are well known from labora- 
tory measurements. T h e same applies to Sn II as a result of the experim ental results of 



Schectman et al.l (120001 ) . IMortonl (120031 ) recommends the theoretical value of iBiemont et al. 



( 119981 ) for the Ge II line at 1237 A. While no more recent theoretical or experimental de- 
terminations are available at the present time, the dichotomy noted above would still be 
present if the results on interstellar Ge were not included. 
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A summary of the depletion patterns for the five elements appears in the top panel of Fig. 
3. Observational results for the amount of depletion onto gra ins both for low and high density 
clouds is sho wn. Each survey used t o construct this figure [Cartledge. Meyer, fc Lauroesch 
J2003h for Kr; ICartledge etaJ (bood ) for Ge; our results for Rb: ISofia et al.l Jl999h for Cd and 
Sn] indicate typical observational scatter of about 0.1 dex in the quoted depletion. Although 
T con d is similar for Ge, Rb, Cd, and Sn, the depletion patterns for Cd and Sn are strikingly 
different than those for Ge and Rb, even when the intrinsic scatter in plots of depletion 
versus T cond are considered. Furthermore, there is less depletion for Cd and Sn than for Kr, 
a noble gas. 

A closer look at the production routes for these n-capture nuclides offers a promising 
solution. This is highlighted in the bottom panel of Fig. 3, where the white bar indicates 
production by low-mass stars and the shaded bar indicates massive stars. The solar system 
abundances for the elements Cd and S n are believed to arise predominantly by the main s 



proce ss involving low- mass AGB stars (IBeer et al.lll992l ; iRaiteri et al.lll993l ; lArlandini et al. 



19991 ). and they have reduced levels of depletion in the ISM. On the other ha nd, elements syn 



thesized by the weak s-process and the r-process that occur in massive stars (jBeer et al.lll992 



Raiteri et al.lll993l : lArlandini et al.lll999l : iThe et al.l 120071 : iHeil et al.ll2008r ) appear to show 
interstellar abundances significantly lower than seen in the solar system. This possible re- 
duction in the contribution of massive stars to interstellar abunda nces for n-capture elements 
contra sts with recent results for lighter elements. For instance, iPrzybilla. Nieva. fc Butler 
(120081 ) provide abundances for several elements in unevo lved, early-type B s tars that are 
i ndistinguishable from those found for the Orion Nebula (jEsteban et al.l I2004J ) and the sun 
(lAsplund et al.ll2005l ). The enhancement of 87 Rb in the Rho Ophiuchus Molecular Cloud 
remains puzzling. Our inferences may not be appropriate if Galactic chemical evolution al- 
tered the mix of s- and r-processed material from 4.5 Gyr ago when the solar system formed 
to the ISM of today. While we are not aware of predictions of the evolution of Rb and its 
isotopes, the weak s-process and r-process from Type II SN might be tied to the evolution 
of O or Mg or Si, while the main s-process to products of AGB stars such as Y or Ba. 

Further efforts in a number of areas will likely improve our understanding. Data on 
interstellar Ga, As, and Pb - other n-capture nuclides - are available, but as of now it is 
more difficult to disentangle the effects of nucleosynthesis from depletion onto grains for 
them. Once this is accomplished, the unexpected results on Rb and its isotopes will likely 
be clarified. A more complete set o f theoretical results, covering all isotopes of interest and 
based on modern stellar codes (e.g-. lThe et al.ll2007l ; IHeil et al.ll2008l ). are needed as well. We 
also see a need for theoretical studies into the evolution of these elements and their isotopes. 
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6. Summary and Conclusions 

The interstellar 85 Rb/ 87 Rb ratio was measured along seven lines of sight, o Per, ( Per, 
AE Aur, HD 147889, x Opli, £ Oph, and 20 Aql. Using K I absorption as a guide, the spectra 
were synthesized to find the velocity, column density, and fe-value of each cloud. Most clouds 
had ratios that agreed, within mutual uncertainties, with the solar system value of 2.59. The 
gas toward HD 147889 had a lower 85 Rb/ 87 Rb ratio (1.03 ± 0.21) and a higher elemental 
ratio (8.6 ± 1.3), suggesting that 87 Rb is enhance d in this direction, consistent with the 



previous results for p Oph A (IFederman et al.l 120041 ) . The anomalous ratios appear to arise 



within the Rho Ophiuchus Molecular Cloud. Furthermore, a comparison of interstellar Rb 
to K shows that the gas in the solar neighborhood appears to be significantly underabundant 
in rubidium, even after considering uncertainties in ionization balance and level of depletion. 
An examination of other species, such as Ge, Kr, Cd, and Sn, seems to suggest nuclides 
synthesized by the n-capture processes occuring in massive stars (weak s and r) have lower 
interstellar abundances than expected. Larger samples for several n-capture elements, from 
both observation and theory, will help clarify the cause for this. 

We thank Dr. Y. Sheffer for his help with his program ISMOD. We acknowledge the use 
of the DIB database of Professor D. York at the University of Chicago, from which promising 
directions were noted. We also thank Prof. P. Boisse for the AE Aur K I A4047 spectrum. 
DLL thanks the Robert A. Welch Foundation of Houston for their support through grant 
number F-634. We made use of the SIMBAD database, operated at Centre de Donees 
Astronomiques de Strasbourg, Strasbourg, France. 
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Fig. 1. — Synthesis of Rb I toward ( Per, x Oph, and HD 147889 based on one-component 
fits. The histogram indicates the data, while the solid line is the fit. The residuals (data 
minus fit) are offset to 1.005, except for HD 147889, where they are offset to 1.01. The 
vertical solid lines above the spectrum of £ Per indicate 85 Rb while the dashed lines indicate 
87 Rb. The lengths show the relative line strengths. The feature in x Oph at 7800.68 A is a 
glitch that did not compromise our spectrum. Note the expanded scale for HD 147889. The 
bottom left panel is fit with a 85 Rb/ 87 Rb ratio of 1.05, while the bottom right panel is fit 
with the solar system ratio; the isotope ratio is given in the lower right corner. 
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Fig. 2. — Synthesis of Rb I toward o Per, AE Aur, ( Oph and 20 Aql based on two-component 
fits. Same as Figure 1. 
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Fig. 3. — The upper panel depicts depletion for warm gas (squares) and cold gas (x's) for 
Ge, Kr, Rb, Cd, and Sn along with their condensation temperatures. The uncertainties for 
the measured abundances are smaller than the symbols. The two groupings in production 
mainly from massive stars and mainly from low-mass stars are ordered by condensation 
temperature. The bottom panel shows production from massive stars (shaded bars) and low 
mass stars (white bars). 



Table 1. Observational Data 



HD 


AT 

Name 


Type 


T T 

V 


a 


Date Observed 


Exposure Time 


T T 11 1 

Unreddened 








(mag) 


(pc) 




(mm) 


Star 


23180 


o Per 


Bl III 


4.98 


450 


Dec. 2005, Jan. 2007 


1470 


a Leo, 7 Cas 


24398 


C Per 


Bl lab 


2.88 


300 


Dec. 2005, Jan. 2007 


946 


a Leo, 7 Cas 


34078 


AE Aur 


09.5 Ve 


6.00 


450 


Dec. 2005, Jan. 2007 


1500 


a Leo, 7 Cas 


147889 




B2 III/IV 


7.95 


135 


Jun. 2006 


560 


a Vir 


148184 


X Oph 


B2 V 


4.42 


150 


May 2005 


810 


a Vir, a Lyr 


149757 


C Oph 


09 V 


2.58 


140 


May 2005, Jun. 2006 


801 


a Vir, a Lyr 


179406 


20 Aql 


B3 V 


5.36 


375 


May 2005, Jun. 2006 


1055 


a Vir, a Lyr 
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Table 2. Rb I Spectra 



Name 


S/N 


W x /a 


o Per 


3300 


8.9 


C Per 


4500 


10.8 


AE Aur 


1500 


10.5 


HD147889 


700 


13.2 


X Oph 


2100 


11.6 


COph 


4100 


16.9 


20 Aql 


2100 


18.4 
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Table 3. Rb I Hyperfine Structure Component Wavelengths and /-values 



Species A (A) /-value 

85 Rb 7800.232 0.290 

7800.294 0.406 

87 Rb 7800.183 0.261 

7800.321 0.435 
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Table 4. Isotope Ratios 



Star 


Isotope 


Vlsr 


N 






b- value 


85 Rb/ 87 Rb 


Reduced 






(km s _1 ) 


(cm- 2 ) 






(km s" 1 ) 


Ratio 


x 2 


o Per 


87 


4.05 


<1.00 x 


10 8 




0.84 








85 


4.05 


(2.72 ± 0.35; 


X 


10 8 


0.84 


>2.72 






87 


7.07 


(2.70 ± 0.30; 


X 


10 8 


1.06 








85 


7.07 


(6.40 ± 0.7i; 


X 


10 8 


1.06 


2.37 ± 0.79 


0.72 


C Per 


87 


5.88 


(3.29 ± 0.29; 


X 


10 8 


1.46 








85 


5.88 


(7.49 ± 0.67; 


X 


10 8 


1.46 


2.28 ± 0.59 


1.75 


AE Aur 


87 


5.34 


<3.00 x 


10 8 




1.03 








85 


5.34 


(8.23 ± 1.07; 


X 


10 8 


1.03 


>2.74 






87 


7.30 


(4.50 ± 0.59; 


X 


10 8 


1.00 








85 


7.30 


(1.13 ± 0.15; 


X 


10 9 


1.00 


2.52 ± 1.03 


1.06 


HD 147889 


87 


2.04 


(4.21 ± 0.35; 


X 


10 9 


1.33 








85 


2.04 


(4.33 ± 0.36; 


X 


10 9 


1.33 


1.03 ± 0.21 


1.44 


X Oph 


87 


0.13 


(4.85 ± 0.38; 


X 


10 8 


0.46 








85 


0.13 


(1.27 ± O.li; 


X 


10 9 


0.46 


2.62 ± 0.63 


2.06 


COph 


87 


-1.21 


(2.99 ± 0.25; 


X 


10 8 


0.45 








85 


-1.21 


(4.72 ± 0.39; 


X 


10 8 


0.45 


1.58 ± 0.32 






87 


-0.17 


(1.54 ± 0.2i; 


X 


10 8 


0.59 








85 


-0.17 


(3.19 ± 0.43; 


X 


10 8 


0.59 


2.07 ± 0.78 


0.94 


20 Aql 


87 


2.06 


(7.26 ± 0.53; 


X 


10 8 


1.28 








85 


2.06 


(2.05 ± 0.15; 


X 


10 9 


1.28 


2.82 ± 0.68 






87 


3.04 


<2.42 x 


10 8 




1.19 








85 


3.04 


(6.06 ± 0.3i; 


X 


10 8 


1.19 


>2.51 


0.90 



-24- 



Table 5. Abundance Ratios from the Neutral Species 



Star 


Vlsr 
(km s _1 ) 


8f 


Rb/ 8 'Rb 


Rb/K a 


85 Rb/K a 


8 'Rb/K a 


o Per 


4.05 




>2.72 


10.0 ± 1.0 


7.3 


± 


1.3 


<2.7 




7.07 


2 


37 ± 0.79 


4.0 ± 1.1 


2.8 


± 


0.7 


1.2 ± 0.4 


C Per 


5.88 


2 


28 ± 0.59 


6.1 ± 0.6 


4.2 


± 


0.3 


1.9 ± 0.3 


AE Aur 


5.34 




>2.74 


3.0 ± 0.6 b 


2.2 


± 


0.8 b 


<0.8 b 




7.30 


2 


52 ± 1.03 


6.1 ± 1.3 b 


4.3 


± 


0.9 b 


1.7 ± 0.4 b 


HD 147889 


2.04 


1 


03 ± 0.21 


8.6 ± 1.3 C 


4.4 


± 


0.7 C 


4.3 ± 0.6 C 


X Oph 


0.13 


2 


62 ± 0.63 


5.4 ± 0.7 


3.9 


± 


0.4 


1.5 ± 0.3 


COph 


-1.21 


1 


58 ± 0.32 


5.7 ± 0.6 


3.5 


± 


0.3 


2.2 ± 0.3 




-0.17 


2 


07 ± 0.78 


6.4 ± 1.1 


4.3 


± 


0.7 


2.1 ± 0.5 


20 Aql 


2.06 


2 


82 ± 0.68 


7.0 ± 0.7 


5.2 


± 


0.5 


1.8 ± 0.3 




3.04 




>2.51 


5.2 ± 0.7 


3.7 


± 


0.8 


<1.5 


p Oph A 




1 


21 ± 0.30 d 


13.0 ± 3.0 d 


7.3 


± 


1.2 d 


6.1 ± 1.3 d 


Solar System 






2.59° 


17.8 ± 3.7° 


12.9° 


4.9 e 



a The values for Rb/K, 85 Rb/K, and 87 Rb/K are in units of 10 
b From data provided by P. Boisse. 



: K I values are from lWelty k, Hobbd (1200 ll ). 



' IFederman et all (120041 ) . 



Loddersl (120031 ). 



